Determination of the Strong Phase in D° —* K + ir Using Quantum-Correlated 

Measurements 
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We exploit the quantum coherence between pair-produced D° and D° in i/>(3770) decays to study 
charm mixing, which is characterized by the parameters x and y, and to make a first determination of 
the relative strong phase 8 between D° — > K + tv~ and D° — » K + tv~ . Using 281 pb _1 of e + e~ collision 
data collected with the CLEO-c detector at E cm — 3.77 GeV, as well as branching fraction input and 
time-integrated measurements of Rm = (x 2 + y 2 )/2 and Rws = T(D° — > K + tc~)/T(D° — » K + tt~) 
from other experiments, we find cos 5 — 1.03^017 ± 0.06, where the uncertainties are statistical and 
systematic, respectively. By further including other mixing parameter measurements, we obtain an 
alternate measurement of cosi5 = 1.10 ± 0.35 ± 0.07, as well as xsin<5 = (4.4l^'g ± 2.9) x 10 -3 and 
5 = (22l^±» 1 )°. 

PACS numbers: 12.15.Ff,13.20.Fc,13.25.Ft,14.40.Lb 



The phenomenon of charm mixing is conventionally 
described by two small parameters, x = (M2 — M\)/Y 
and y = (T2 — Pi)/2r, where M\ t 2 and Ti.2 are the 
masses and widths, respectively, of the CP-odd (Di) 



and CP-even (D2) neutral D meson mass eigenstates, 
and r = (ri + T2)/2. Many previous searches for charm 
mixing have used D° decay times to attain first-order 
sensitivity to y. Lifetimes of D° decays to CP eigen- 
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states determine y, while doubly Cabibbo-suppressed 
(DCS) transitions probe Rm = {x 2 + y 2 )/2 and a mode- 
dependent quantity, y'. For the most widely used DCS 
mode, D° — ► K +/ n~ , y' = ycosS — a;sin<5, where —S is 
the relative phase between the DCS amplitude and the 
corresponding Cabibbo-favored P/° — ► K + ?r~ amplitude: 
(K+7r-\D°)/(K+Tr-\D ) = re~ lS . We adopt a conven- 
tion in which 8 corresponds to a strong phase, which 
vanishes in the SU(3) limit To date, 8 has not been 
measured, so measurements of y and y' have not been 
directly comparable. The magnitude r of the amplitude 
ratio is approximately 0.06. 

In this Letter, we implement the method described in 
Ref. for measuring y and cos 8 using quantum correla- 
tions at the V(3770) resonance [H, where D°D° pairs 
produced in e + e~ collisions are in a C-odd eigenstatc. 
We extract these parameters from decay rates to single 
tags (ST), which are individually reconstructed D° or 
D° candidates, and double tags (DT), which are events 
where both D° and D° are reconstructed. CP violation 
in D and K decays are negligible second order effects. 

To first order in x and y, the rate T DO j}o(i, j) for C- 
odd D°D° decay to final state {i,j} follows from the 
anti-symmetric amplitude M.if 

T DOB o(i,j) oc M% = \AiAj-AiAjl 2 

= |(i|D 2 )(i|A)-(i|A)0-|i) 2 )| 2 , (i) 

where Aj = (i\D°), A\ = (i\D°), and we have used 
|£>i) = [\D°) ± \D°)]/V2. Using S± and e ± to de- 
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note CP± eigenstates and semileptonic final states, re- 
spectively, these amplitudes are normalized such that 
B K -„+ « A%-_ n+ (l + ry cos 8 + rxwa.5), B s± » A| ± (1=f 
y), and S e ~ A 2 . Quantum correlations affect neither 
the total D°D° rate (and hence the number AT of P°P/° 
pairs produced) nor the ST rates. DT final states with 
pairs of CP eigenstates, however, are affected maximally; 
same-CP {S±, S±} states are forbidden, while opposite- 
CP {S+,5_} states are doubled in rate relative to un- 
correlated decay. In general, the correlations introduce 
interference terms that can depend on y and 8. 

D°D° decay involving a final CP eigenstate naturally 
selects the D\D^ basis. As a result, the branching frac- 
tion for an associated semileptonic decay probes y. While 
the semileptonic decay width itself does not depend on 
the CP eigenvalue, the total width of the parent D\ or 
Di meson does: Ti = T(l =p y). Thus, the D\ semilep- 
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tonic branching fraction is B e /(1 =F y), and the effective 
quantum-correlated D°D° branching fraction (J rcor ) for 
a {S±,e} final state is Fg™ e ~ 2Bs ± B e (l±y), where the 
factor of 2 arises from the sum of e + and e~ rates. When 
combined with estimates of B e and Bs ± from ST yields, 
external sources, and flavor-tagged semileptonic yields, 
this equation allows y to be determined. 

If an S+ and a K~tt + decay occur in the same event, 



then the K ir + was produced by a Pi, and .F|° r K „ is 

FsIk, = |<S + |P 2 )<^ + |Pi)| 2 
= A 2 s+ \A k -k+ + A K -^+\ 2 
= A 2 + A 2 K ^ + \l + re-* s \ 2 
w B S+ B K - v+ (1 + R ws + 2r cos S + y), (2) 

where Pws is the wrong-sign rate ratio, which depends 
on x and y because of the interference between DCS and 
mixing transitions: P ws = r(P° -> K-tt+)/T{D [) -> 
K-TT+) = r 2 + ry' + R M . Similarly, the {S-,Kn} DT 
yield probes Bs_ B K -^+ (1 + Pws — 2r cos 5 — y), and the 
asymmetry between these two DT yields gives cos 8, given 
knowledge of Bs ± , r, and y. 

Table U shows J rcor for all categories of final states con- 
sidered in this analysis: , S±, and e ± . Comparison 
of jF cor with the uncorrelated effective branching frac- 
tions, T unc , also given in Table HI provides rcosS, y, r 2 , 
x , and rxsinS. These five parameters are extracted by 
combining our ST and DT yields with external branch- 
ing fraction measurements in a least-squares fit [1]. The 
external measurements, from incoherently produced D° 
mesons, provide one measure of Bi. The ST event yields 
provide a second measure; since each event has one D° 
and one D , inclusive rates correspond to uncorrelated 
branching fractions. The fit averages these estimates, 
and we extract updated Bi. Finally, the DT/ST com- 
parison provides Af, so the fit requires no knowledge of 
luminosity or D°D° production cross sections. 

TABLE I: Correlated (C-odd) and uncorrelated effective 
D°D° branching fractions, T cor and T unc , to leading order 
in x, y, and Rws, divided by Bi for ST modes i (first sec- 
tion) and BiBj for DT modes {i,j} (second section). Charge 
conjugate modes are implied. 



Mode 


Correlated 


Uncorr. 


K~n + 


1 + Rws 


f + i?WS 


s± 


2 


2 


K-TT + ,K-TT+ 


Rm 


-Rws 




(f + Rws) 2 — 4r cos<5(r cos 5 - 


h y) 1 + R W s 


K-TT+,S± 


f + -Rws ± 2r cos 8 ± y 


f + Rws 


K~Tv + ,e- 


1 — r y cos 8 — rx sin 8 


f 


S±,S± 





1 


S+,S- 


4 


2 


S±,e~ 


l±y 


f 



We analyze 281 pb _1 of e + e~ collision data produced 
by the Cornell Electron Storage Ring (CESR) at E cm = 
3.77 GeV and collected with the CLEO-c detector, which 
is described in detail elsewhere [6j. We reconstruct the 
D° and D° final states listed in TableHU with tt% -> 77, 
lo — * 7T + 7r~7r°, and Kg — > 7r + 7r~. Signal and background 
efficiencies, as well as probabilities for misreconstruct- 
ing a produced signal decay in a different signal mode 
(crossfeed), are determined from simulated events that 
are processed in a fashion identical to data. 
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TABLE II: D final states reconstructed in this analysis. 



Type 


Final States 


Flavored 


K~n+, K+tr- 


s+ 


K+K~, 7T + 7V-, K° S 7Y°7V , K°n° 


S- 




e± 


Inclusive Xe + v e , Xe 



TABLE III: ST and DT yields, efficiencies, and their statis- 
tical uncertainties. For DT yields, we sum groups of modes 
and provide an average efficiency for each group; the number 
of modes in each group is given in parentheses. Modes with 
asterisks are not included in the standard and extended fits. 



Mode 






Yield 


Efficiency (%) 








25374 ± 168 


64.70 ± 0.04 


K+tt- 






25842 ± 169 


65.62 ± 0.04 


K+K' 






4740 ± 71 


57.25 ± 0.09 


7T + 7T~ 






2098 ± 60 


72.92 ± 0.13 


K g-K 7T 




2435 ± 74 


12.50 ±0.06 








7523 ± 93 


29.73 ± 0.05 


K° S V 






1051 ± 43 


10.34 ±0.06 


K° s uj 






3239 ± 63 


12.48 ± 0.04 




,K+w* 


(2) 


4±2 


40.2 ± 2.4 


K--K+ 




(1) 


600 ± 25 


41.1 ± 0.2 




,S+ (8) 




605 ± 25 


26.1 ± 0.1 




,S- (6) 




243 ± 16 


12.3 ± 0.1 




,e T (2) 
(9*) 




2346 ± 65 


45.6 ± 0.1 


s + ,s + 




10 ±6 


12.5 ± 0.6 


S-,S- 


(6*) 




2 ± 2 


3.9 ±0.2 


S+,S- 


(12) 




242 ± 16 


7.7 ±0.1 


S+,e* 


(6) 




406 ± 44 


22.2 ± 0.1 


S-,e* 


(6) 




538 ± 40 


13.8 ±0.1 



The D candidate selection and yield determination 
procedures are described in a companion article 0] and 
are summarized below. Hadronic final states without 
K\ mesons are fully reconstructed via two kinematic 
variables: the beam-constrained candidate mass, M = 
^/Eq/c 4 — p D I c 2 , where is the D° candidate mo- 
mentum and Eo is the beam energy, and AE = Ed — Eq, 
where Ed is the sum of the D° candidate daughter en- 
ergies. We extract ST and DT yields from M distribu- 
tions using unbinned maximum likelihood fits (ST) or by 
counting candidates in signal and sideband regions (DT). 

Because most H\ mesons and neutrinos produced at 
CLEO-c are not detected, we only reconstruct modes 
with these particles in DTs, where the other D in the 
event is fully reconstructed. Ref. Q describes the miss- 
ing mass technique used to identify K^ir candidates. 
For semileptonic decays, we use inclusive, partial recon- 
struction to maximize efficiency, demanding only that 
the electron be identified. Electron identification utilizes 
a multivariate discriminant Q that combines measure- 
ments from the tracking chambers, the electromagnetic 
calorimeter, and the ring imaging Cerenkov counter. 

Table IIIII gives yields and efficiencies for 8 ST modes 
and 58 DT modes, where the DT modes have been 



grouped into categories. Fifteen of the DT modes are for- 
bidden by CP conservation and are not included in the 
nominal fits. In general, crossfeed among signal modes 
and backgrounds from other D decays are smaller than 
1%. Modes with KgTr°ir° have approximately 3% back- 
ground, and yields for {/ < iT =F 7r ± , X T 7r ± } and {S±, S±} are 
consistent with being entirely from background. 

External inputs to the standard fit include measure- 
ments of i?M, -Rws, Bk-k+i an d £>s±j as well as an in- 
dependent B K o^o from CLEO-c, as shown in Table H*Vl 
i?ws is required to constrain r 2 , and thus, to convert 
r cos 8 and rx sin i5 to cos 5 and x sin 5. We also perform 
an extended fit that uses the external mixing parameter 
measurements shown in Table [V] These fits incorporate 
the full covariance matrix for these inputs, accounting for 
statistical overlap with the yields in this analysis. Covari- 
ance matrices for the fits in Ref. [HI have been provided 
by the CLEO, Belle, and BABAR collaborations. 



TABLE IV: Averages of external measurements used in the 
standard and extended fits. Charge-averaged D° branching 
fractions are denoted by final state. 



Parameter 



Average 



-Rws 

Rm 

R-TI+ 

K~K+ 

ir~ tv + I K~ tv + 

K L TT 



0.00409 ± 0.00022 [10] 
0.00017 ± 0.00039 [11] 
0.0381 ± 0.0009 [12] 
0.1010 ±0.0016 [13] 
0.0359 ± 0.0005 [13] 
0.0097 ± 0.0003 [8] 
0.0115 ±0.0012 [12] 
0.00380 ± 0.00060 [12] 
0.0130 ± 0.0030 [12] 



TABLE V: Averages of external measurements used only in 
the extended fit. 



Parameter 


Average 


y 


0.00662 ± 0.00211 [13, 14, 15] 


X 


0.00811 ±0.00334 [15] _ 


r 2 


0.00339 ± 0.00012 [16] 


y' 


0.0034 ± 0.0030 [16] 


x' 2 


0.00006 ± 0.00018 [16] 



Systematic uncertainties include those associated with 
efficiencies for reconstructing tracks, Kg decays, ir° de- 
cays, and for hadron identification (see Refs. [1, 0]). 
Other sources of efficiency uncertainty include: AE re- 
quirements (0.5-5.5%), T) reconstruction (4.0%), electron 
identification (1.0%), modeling of particle multiplicity 
and detector noise (0.1-1.3%), simulation of initial and fi- 
nal state radiation (0.5-1.2%), and modeling of resonant 
substructure in K s tt°tt° (0.7%). We also include addi- 
tive uncertainties of 0.0-0.9% to account for variations 
of yields with fit function. 
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These systematic uncertainties are included in the co- 
variance matrix given to the fitter, which propagates 
them to the fit parameters. The other fit inputs deter- 
mined in this analysis are ST and DT yields and efficien- 
cies, crossfeed probabilities, background branching frac- 
tions and efficiencies, and statistical uncertainties on all 
of these measurements. Quantum correlations between 
signal and background modes are accounted for using as- 
sumed values of amplitude ratios and strong phases that 
are systematically varied and found to have negligible ef- 
fect. Using a simulated C-odd D°D° sample 15 times 
the size of our data sample, we validated our analysis 
technique by reproducing the input branching fractions 
and mixing parameters with biases due to our procedures 
that were less than one-half of the statistical errors on the 
data and consistent with zero. 

Table |VT] shows the results of the data fits, excluding 
the 15 same-CP DT modes. Our standard fit includes 
the measurements in Table IIVI but not Table [V] In this 
fit, xsin<5 is not determined reliably, so we fix it to zero, 
and the associated systematic uncertainty is ±0.03 for 
cos 8 and negligible for all other parameters. We obtain 
a first measurement of cos<5, consistent with being at the 
boundary of the physical region. Our branching fraction 
results do not supersede other CLEO-c measurements. 



TABLE VI: Results from the standard fit (with Table |TV] in- 
puts) and the extended fit (with Table ITV} fVl inputs) . Uncer- 
tainties are statistical and systematic, respectively. Charge- 
averaged D° branching fractions are denoted by final state. 



Parameter 



Standard Fit 



Extended Fit 



N (10°) 
v (io- 3 ) 

r 2 (1(U 3 ) 
cos S 

x 2 (IO" 3 ) 
x sin (5 (1(U 3 ) 
K--K+ (%) 
K~K+ (1(U 3 ) 
TY-n+ (10 -3 ) 
K° s n°iT (1(U 3 ) 

K°s*° {%) 
K° S V (10- 3 ) 
K° s co (%) 
X~e + u e (%) 
K° L n° (%) 
Xflt/ndof 



1.042 ±0.021 ±0.010 
-45 ± 59 ± 15 
8.0 ± 6.8 ± 1.9 
1.03 ±0.19 ±0.06 
-1.5 ±3.6 ±4.2 

(fixed) 
3.78 ±0.05 ±0.05 
3.87 ±0.06 ±0.06 
1.36 ±0.02 ±0.03 
8.34 ±0.45 ±0.42 
1.14 ±0.03 ±0.03 
4.42 ±0.15 ±0.28 
1.12 ±0.04 ±0.05 
6.54 ±0.17 ±0.17 
1.01 ±0.03 ±0.02 



1.042 ±0.021 ±0.010 
6.5 ±0.2 ±2.1 
3.44 ±0.01 ±0.09 
1.10 ±0.35 ±0.07 
0.06 ±0.01 ±0.05 

4.4 ±2.4 ±2.9 
3.78 ±0.05 ±0.05 
3.88 ±0.06 ±0.06 
1.36 ±0.02 ±0.03 
8.35 ±0.44 ±0.42 
1.14 ±0.03 ±0.03 
4.42 ±0.15 ±0.28 
1.12 ±0.04 ±0.05 
6.59 ±0.16 ±0.16 
1.01 ±0.03 ±0.02 



30.1/46 



55.3/57 



The likelihood curve for cos 8, shown in Fig.[T^,, is com- 



puted as C 



-(x 2 



at various fixed values of 



cos 8. It is highly non-Gaussian, so we assign asymmet- 
ric uncertainties (which still do not fully capture the non- 
linearity) by finding the values of cos<5 where A\ 2 = 1 
to obtain cosiJ = 1.03^°'^ ± 0.06. This non-linearity 
stems from the use of r cos 8 to determine cos 8, which 
causes the uncertainty on cos 8 to scale roughly like 1/r. 
Because r 2 is obtained from i?wsj an upward shift in y 



lowers the derived value of r 2 (for positive rcosS), and 
the resultant uncertainty on cos 8 increases, as illustrated 
by Fig. [lb. For values of |cos<5| < 1, we also compute C 
as a function of \8\, and we integrate these curves within 
the physical region to obtain 95% confidence level (CL) 
limits of cos<5 > 0.07 and \8\ < 75°. 



3950108-001 




cos5 

FIG. 1: Standard fit likelihood including both statistical and 
systematic uncertainties for cos 8 (a), and simultaneous like- 
lihood for cos 5 and y (b) shown as contours in increments of 
la, where a = yj Ax 2 - The hatched region contains 95% of 
the area in the physical region. 



3950108-002 




40 

8 (degrees) 



FIG. 2: Extended fit likelihood including both statistical and 
systematic uncertainties for cos 8 (a), x sin 8 (b), 8 (c), and 
simultaneous likelihood for cos 8 and x sin 8 (d) shown as con- 
tours in increments of la, where a = y/ Ax 2 - The hatched 
regions contain 95% of the area in the physical regions. For 
8, the fit fails to converge beyond the limits of the plot. 

When combined with previous measurements of y and 
y' , our measurement of cos 8 also gives xsmS. Table IVT1 
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shows the results of such an extended fit that includes 
external inputs from both Table IIVI and Table [V] The 
resultant value of y includes the CLEO-c measurement 
from the standard fit, but the precision is dominated 
by the external y measurements. The overall uncer- 
tainty on cos 5 increases to ±0.36 because of the non- 
linearity discussed above. However, unlike the standard 
fit, the likelihood for cos<5 is nearly Gaussian, as shown 
in Fig. [2k- The correlation coefficient between cos 6 and 
xsin(5 is 0.56, and we assign asymmetric uncertainties of 
xainS = {4A±f 7 8 ± 2.9) x 10 -3 . By repeating the fit at 
various simultaneously fixed values of cos 5 and sin 5, we 
also determine 6 — (22^^-11 )°- The corresponding 95% 
CL intervals within the physical region are cos<5 > 0.39, 
x sin S G [0.002,0.014], and 5 G [-7°, +61°]. Perform- 
ing this extended fit with y, x 2 , and xsinJ fixed to zero 
results in a change in x 2 of 25.1, or a significance of 5. Oct. 

By observing the change in l/tr^ as each fit input is 
removed, we identify the major contributors of informa- 
tion on y to be the {S±,e} yields (90%) and {Kir, e} 
yields (10%). For cos<5, the {Kir, S±\ DT yields and the 
ST yields simultaneously account for 100%. We also find 
that no single input or group of inputs exerts a pull larger 
than 3cr on cos 8 or y. Moreover, removing all external 
inputs gives branching fractions consistent with those in 
Table IIVI Finally, if we determine y only from K + K~ 
and 7r + 7r~ input, as in previous direct measurements, the 
result is consistent with the value in Table [VTl 

We also allow for a C-even D°D° admixture in the 
initial state, which is expected to be 0(1O -8 ) by 
including the 15 {S±,S±} DT yields in the fit. These 
modes limit the C-even component, which can modify 
the other yields as described in Ref. Q. In both the 
standard and extended fits, we find a C-even fraction 
consistent with zero with an uncertainty of 2.4%, and 
neither the fitted parameters nor their uncertainties are 
shifted noticeably from the values in Table PVTl 

In summary, using 281 pb _1 of e + e~ collisions pro- 
duced at the ip(3770), we make a first determination of 
the strong phase 5, with cos<5 = 1.03^q]i7 ± 0.06. By 
further including external mixing parameter measure- 
ments in our analysis, we obtain an alternate measure- 
ment of cos 5 = 1.10 ± 0.35 ± 0.07, as well as x sin 6 = 
(4.4^-7 ± 2.9) x 10~ 3 and S = (22±\l± 9 n )°. 
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